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The X-ray spectrum of the nearest ultraluminous X-ray source, M33 X-8, obtained by Suzaku during 
\ 2010 January 11 - 13, was closely analyzed to examine its nature. It is, by far, the only data with the highest 

signal statistic in 0.4 - 10 keV range. Despite being able to reproduce the X-ray spectrum, Comptonization 
of the disk photons failed to give a physically meaningful solution. A modified version of the multi-color 
disk model, in which the dependence of the disk temperature on the radius is described as r~ p with p being 
a free parameter, can also approximate the spectrum. From this model, the innermost disk temperature 



and bolometric luminosity were obtained as Tj n = 2.001qq5 keV and Ldisk = 1-36 x 10 39 (cosi) 1 ergs s 1 , 
<*S { respectively, where i is the disk inclination. A small temperature gradient of p = 0.535^'^, together with 

the high disk temperature, is regarded as the signatures of the slim accretion disk model, suggesting that 
M33 X-8 was accreting at high mass accretion rate. With a correction factor for the slim disk taken into 



account, the innermost disk radius, R[ n = 81.9^6 5( cos *) °" 5 km, corresponds to the black hole mass of 



M ~ lOM (cosi) -0 - 5 . Accordingly, the bolometric disk luminosity is estimated to be about 80(cosi) _a5 % 







of the Eddington limit. A numerically calculated slim disk spectrum was found to reach a similar result. 
Thus, the extremely supcr-Eddington luminosity is not required to explain the nature of M33 X-8. This 
conclusion is utilized to argue for the existence of intermediate mass black holes with M > 100A1© radiating 
at the sub/trans-Eddington luminosity, among ultraluminous X-ray sources with Ldisk ^ 10 40 ergs s . 

Introduction lar mass. In the second alternative, a super Eddington 

luminosity is permitted. Consequently, the ULXs are re- 
A number of nearby galaxies are known to host hum- garded as a stellar-mass black hole (M<10M Q ), accreting 



in 

^ ! nous X-ray sources with an X-ray luminosity of L x > 10 39 a t a super-critical mass accretion rate (M > L E /c 2 with 

ergs s _1 (e.g., Fabbiano & Trinchieri 1987; Swartz et al. c being the speed of light). 

2011). They are usually referred to as ultraluminous X- i t [ s f vital importance to construct a reliable tech- 
ray sources (ULXs; Makishima et al. 2000). The ULXs mque to gauge the mass and/or mass accretion rate of the 
^ ; are found to be less common in the Local Universe, as no ULXs, in order to distinguish these two conflicting inter- 
persistent ULX is yet identified in our Galaxy and M31. pre tations. For the purpose, the dynamical mass determi- 
In spite of significant efforts in more than three decade na tion is thought to provide inevitably the most stringent 
since their discovery, we have not yet arrived at the gen- constraint, as is the case of the Galactic black hole binaries 
eral conclusion on their nature. (e.g., Shahbaz et al. 1999; Hynes et al. 2003). However, 
It is widely believed from their observational features the dynamical mass of the ULXs has remained yet to be 
(e.g., Makishima et al. 2000) that the ULX phenomena explored. 

are related to accreting black holes. Because the lumi- X-ray spectra and their variability have been one of 

nosity of an accreting black hole is though to increase as the most valuable tools to probe the properties of the 

a function of the black hole mass M and accretion rate ULXs. A naive inspection on the X-ray spectra in the 

M on the simplest assumption, two distinct ideas are cur- q.5 - 10 keV range categorizes the ULXs into two types; 

rently proposed to interpret the high X-ray luminosity of those with an X-ray spectrum dominated by emission from 

the ULXs. The first option supposes that the luminosity accretion disk which is approximated by the multi-color 

of the ULXs should not exceed the Eddington luminosity disk (MCD; Mitsuda et al. 1984; Makishima et al. 1986) 

(£ E ), as most Galactic black hole binaries meet this cri- mo del, and those exhibiting a power-law (PL) like X-ray 

terion (e.g., Kubota & Makishima 2004). As a result, the spectrum. In addition, a growing number of ULXs are 

ULXs arc required to harbor an intermediate-mass black rC p rted to make a spectral transition between the MCD- 
holes with a mass of M ^> 10A/© with M© being the so- 
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like and PL-like states (e.g., Kubota et al. 2001b; Isobe 
et al. 2009). These spectral behaviors apparently resem- 
ble those of Galactic black hole binaries in the classical 
high/soft and low/hard states. However, several observa- 
tional obstacles are known for this assumption as follows 
(e.g., Makishima et al. 2000; Mizuno et al. 2001). The 
X-ray spectrum of the MCD-like ULXs requires a flat- 
ter radial temperature distribution within the disk (e.g., 
Tsunoda et al. 2006) in comparison with that in the stan- 
dard accretion disk (Shakura & Sunyaev 1973) which is 
adopted in the MCD model. The innermost disk ra- 
dius i?j n of the MCD-like ULX appears to be variable 
as R- m oc Xj^ 1 with Ti n being the innermost disk temper- 
ature (e.g., Mizuno et al. 2001; Isobe et al. 2009), while 
the Galactic black hole binaries in the high/soft state are 
reported to follow a constancy of R[ n (e.g., Ebisawa et al. 
1993) which indicates that i? in corresponds to the last sta- 
ble circular orbit of the accretion disk. The photon index 
of the PL-like ULXs (r > 2; e.g., Kubota et al. 2002; Isobe 
et al. 2009) tends to be larger than those of the Galactic 
black hole binaries in the low/hard state (r ~ 1.5; Tanaka 
1997). 

Recent progresses in observational (e.g., Kubota & 
Makishima 2004; Kubota & Done 2004) and theoretical 
(e.g., Watarai et al. 2000) studies on Galactic black hole 
binaries and their accretion disk, especially near the high- 
est end of their X-ray luminosity, newly identified two 
spectral states called the very high and slim disk ones, in 
addition to the classical low/hard and high/soft states. In 
the very high state, the Galactic black hole binaries are 
reported to exhibit a steep PL spectrum with a photon 
index of L > 2.4, which is softer than the PL component 
observed in the high/soft state (r = 2 -2.2). The steep PL 
spectrum in this situation is ascribed to Comptonization 
of the disk photons within the hot corona surrounding the 
accretion disk. In contrast, when the Galactic black hole 
binaries become brighter than in the very high state, they 
are reported to enter the slim disk state in which a disk 
dominant X-ray spectrum with a flat disk temperature 
gradient was observed, in addition to an apparent anti- 
correlation between the innermost disk radius and tem- 
perature (Kubota & Makishima 2004). These properties 
are successfully explained by the theory of a slim accre- 
tion disk, where an optically-thick advection and photon 
trapping (Ohsuga et al. 2005) effectively work due to a 
high accretion rate (M > Le/c 2 ). 

Triggered by these findings, it is proposed that the 
MCD-like and PL-like ULXs reside in the slim disk 
and very high states, respectively (e.g., Mizuno et al. 
2001; Kubota et al. 2002; Tsunoda et al. 2006; Mizuno 
et al. 2007; Isobe et al. 2009). If this idea is system- 
atically proved to be convincing, we will have acquired 
an empirical method to restrict the black hole mass of 
the ULXs. This is because the X-ray luminosity of the 
Galactic black hole binaries in these states are reported 
typically as Lx = (0.3 - 1)Le (e.g., Kubota & Makishima 
2004; Kubota et al. 2001a). In other word, this interpre- 
tation consequently indicates that the ULXs are shining 
at a sub- or trans-Eddington luminosity. Independently, 




Fig. 1. Background-inclusive Suzaku XIS image in the range 
of 0.5 - 10 keV of M33 (red) with a bin size of 10" X 10", on 
which the optical image taken with the Digitized Sky Survey 
(blue) is superposed. The data from all the active XIS chips 
(XIS 0, 1 and 3; Koyama et al. 2007) were summed up. The 
XIS image was smoothed with a two-dimensional Gaussian 
kernel of a 20" radius, while it is not corrected for the ex- 
posure. The regions irradiated by the radioactive calibration 
sources at the corners of the CCD chips were removed. The 
X-ray sources detected with ROSAT are indicated with the 
green crosses, accompanied by the source name (Trinchicri et 
al. 1988). The XIS signals from M33 X-8 and background are 
integrated within the solid and dashed circles, respectively. 

based on numerical simulation, Vierdayanti et al. (2008) 
extended to the slim accretion disk the mass estimation 
from the innermost disk radius, which is widely accepted 
for Galactic black hole binaries in the high/soft state ac- 
companied with the standard accretion disk. It is impor- 
tant to confirm observationally the consistency between 
these techniques. 

Recently, it is verified that Suzaku (Mitsuda et al. 2007) 
is very useful to investigate the X-ray spectral proper- 
ties of the ULXs (Mizuno et al. 2007; Miyawaki et al. 
2009; Isobe et al. 2008; Isobe et al. 2009). These results 
motivated us to perform a Suzaku observation of the near- 
est ULX, M33 X-8 (Trinchieri et al. 1988), in which its X- 
ray spectrum was obtained with the highest signal statis- 
tics in the 0.4 - 10 keV range. In the present paper, we 
utilized extensively the the Suzaku data of the source to 
investigate the nature of the ULXs. 

2. Observation and Data Reduction 

Within ~ 20 arcmin from the center of the nearby nor- 
mal galaxy, M33, more than 10 bright X-ray sources were 
discovered with the Einstein observatory (Trinchieri et al. 
1988). Among these sources, X-8 located near the center 
of the galaxy is known to show the ULX characteristics 
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(Makishima et al. 2000). Thanks to it proximity at a 
distance of 795 kpc (van den Bergh 1991), it extensively 
studied with the previous and active instruments includ- 
ing ASCA (Takano et al. 1994), BeppoSAX (Parmar et 
al. 2001), Chandra (La Parola et al. 2003), and XMM- 
Newton (e.g., Weng et al. 2009). The typical X-ray lumi- 
nosity of M33 X-8, Lx ^ 10 39 ergs s _1 , is at the lowest 
end of those of ULXs; only the Galactic microquasar, GRS 
1915+105, is confirmed to reach occasionally such a high 
X-ray luminosity value (e.g., Ueda et al. 2010). Some au- 
thors claimed a flat disk temperature distribution to sug- 
gest the slim disk interpretation for the source (Foschini 
et al. 2006; Weng et al. 2009), although several conflicting 
ideas were proposed (Foschini et al. 2004; Middleton et al. 
2011). These make M33 X-8 one of the most important 
object to disentangle the nature of the ULXs by connect- 
ing the properties of the ULXs with those of the Galactic 
black hole binaries. 

The Suzaku observation of M33 X-8 was performed on 
2010 January 11 - 13. The X-ray imaging spectrome- 
ter (XIS; Koyama et al. 2007) and hard X-ray detector 
(HXD; Takahashi et al. 2007) was operated in the normal 
clocking mode with no window option and in the normal 
mode, respectively The source was placed at the HXD 
nominal position of the X-ray telescope (Scrlcmitsos et al. 
2007) above the XIS. Unfortunately, due to a significant 
increase in the thermal noise of the HXD PIN, the HXD 
effective exposure was found to reduce to only about 6.1 
ks. Therefore, we concentrate on the XIS data in the 
present paper. 

The data reduction was performed with the standard 
software package HEADAS 6.10. We reprocessed the XIS 
data with the HEADAS tool aepipeline, by referring to 
the calibration database (CALDB) as of 2011 February 
10. In accordance with the standard manner, the follow- 
ing filtering criteria were adopted; the spacecraft is outside 
the south Atlantic anomaly (SAA), the time after an exit 
from the SAA is larger than 436 s, the geometric cut-off 
rigidity is higher than 6 GV, the source elevation above 
the rim of bright and night Earth is higher than 20° and 
5° , respectively, and the XIS data are free from telemetry 
saturation. As a result, we derived 91.0 ks of good expo- 
sure. The XIS events with a grade of 0, 2, 3, 4, or 6 were 
picked up for the scientific analysis below. 

3. Results 

3.1. X-ray image 

Figure 1 displays the 0.5 - 10 kcV XIS image of M33 
X-8, which is superposed on the optical image of the 
host galaxy taken from the digitized sky survey (Lasker 
et al. 1990). The brightest X-ray source in the im- 
age, detected at the center of the host galaxy, corre- 
sponds M33 X-8. A systematic offset of 15.6" was no- 
ticed in the XIS coordinate determination, while it is 
found to be within the current systematic uncertainties 
(Uchiyama et al. 2008). We thus re-calibrated the XIS 
coordinate by ourselves, referring to the ROSAT posi- 
tion of M33 X-8 itself ((a 20 oo, £2000) = (01 h 33 m 50f58, 



BGD-subtracted XIS FI lightcurve ofM33 X-8 
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Fig. 2. Background-subtracted XIS FI lightcurve of M33 
X-8. Each time bin is set to 5760 s, corresponding to the or- 
bital period of Suzaku. The top and middle panels show the 
count rate in the soft (0.5 - 2 keV) and hard bands, respec- 
tively. The source hardness, simply evaluated as the hard-to- 
soft count-rate ratio, is shown in the bottom panel. Assuming 
the Galactic absorption toward M33, iVn = 1.1 x 10 21 cm -2 
(Kalbcrla et al. 2005), a simple PL model with a photon index 
of r = 2.0 and 1.8 corresponds to the hardness shown with the 
dotted and dashed lines, respectively. 

+30 d 39 m 36?7); Trinchieri et al. 1988). Within the XIS 
field of view, another X-ray source M33 X-7, a candi- 
date of a massive spinning black hole with a mass of 
M = (15.6±3.2)AY© (Abubekerov et al. 2010), is clearly 
seen at a position consistent with the ROSAT one ((0:2000- 
£2000) = (01 h 33 m 34?08, +30 d 32 m 13?2); Trinchieri et al. 
1988). However, due to the insufficient signal statistics 
from the XIS, we regard its X-ray spectral properties as 
irrelevant to the scope of the present paper. 

In the scientific analysis below, the XIS signals from 
M33 X-8 were accumulated within a solid circle in figure 
1 with a 3' radius centered on the source, while the back- 
ground level was evaluated from a source free region with 
the same radius, denoted by the dashed circle in figure 
1. Although there is a faint ROSAT X-ray source, M33 
X-15, near the edge of the source integration region, we 
estimated its contamination as negligible in comparison 
with the X-ray flux from M33 X-8 (less than 1%), based 
on the ROSAT result (Trinchieri et al. 1988). 

3.2. X-ray lightcurve 

Figure 2 shows the background-subtracted lightcurves 
of M33 X-8 in the soft (0.5 - 2 keV) and hard (2 - 10 keV) 
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Table 1. Summary of the model fitting to the XIS spectra of M33 X-8 



[Vol. , 



Model 


MOD" 


PL" 


MCD+PL 


MCD+THCOMP extended MCD 


Fx* 


1.50 


1.60 


1.57 


1.55 


1.56 




0.57 


0.83 


0.75 


0.65 


0.71 


A H (10 21 cm- 2 ) 


0.0 


2.38 


1.80 ±0.13 




i 47+o.os 
i>4: '-Q.Q7 


H„ (kcV) 


1.27 




i oq+0.07 
i-OO -0.06 


0.36jig 




i?in X 


51.4 




29.2±|f 


<88.0 


8i.9i§ : § 


£disk § 


0.62 




0.29 


<0.01 


1.36 


V 










u.ooo_ 005 


r 




2.17 


o on+0.08 

^• ou -0.07 


1 qi+004 




T c (keV) 








1 62+ 010 




x7dof 


3032.7/1373 2429.7/1373 1535.8/1371 


1549.6/1370 


1541.8/1372 



* Absorption-inclusive 0.5 - 10 keV flux in 10 ergs cm t 
t Absorption-corrected 0.5 - 10 keV luminosity in 10 39 (cosi)~ 
t The true innermost disk radius in (cosi) -0 - 5 km. The spectral hardening factor of ft = 1.7 and the correction 
factor for the boundary condition of £ = 0.412 are applied in the case of the MCD component, while the extended 
MCD model adopted n = 3 (Watarai & Mineshige 2003) and £ = 0.353 (Vierdayanti et al. 2008). 



Absorption-corrected bolometric disk luminosity in 10 39 (cosi) 
The errors are omitted, because the fit was unacceptable. 



ergs s 



■\ evaluated in 0.001 - 100 keV 



bands, where the data from the two front-side illuminated 
(FI) CCD chips (XIS and 3; Koyama et al. 2007) were 
summed up. The source hardness, plotted in the bottom 
panel of figure 2, was simply calculated as the ratio of the 
hard band count rate to the soft band one. 

Summed over the two FI CCDs, we measured the time- 
averaged signal count rate of M33 X-8 as 0.77 ± 0.03 
cts s^ 1 and 0.56 ±0.03 cts s _1 , respectively, in the soft 
and hard band. Figure 2 suggests a weak but statis- 
tically significant intensity variation (x 2 /dof = 145.3/37 
and 142.9/37) during the Suzaku exposure in the soft and 
hard ranges with a standard deviation of 4.8% and 5.4%, 
respectively, around the average count rate. The hardness 
of the object was revealed to be relatively unchanged dur- 
ing the Suzaku exposure, although the source appears to 
exhibit a possible spectral softening near the end of the 
exposure. 

3.3. X-ray spectrum 

Since the change of hardness of M33 X-8 was suggested 
to be insignificant during this observation as shown in fig- 
ure 2, we consider only the time- averaged spectrum of the 
source. Figure 3 shows the background-subtracted XIS 
spectrum of M33 X-8 in the range of 0.4 - 10 keV, derived 
with the FI and back-side illuminated (BI; XIS 1, Koyama 
et al. 2007) CCD chips, without removing the instrumen- 
tal response. Both the FI and BI spectra were binned into 
energy intervals each with at least 100 events. In compar- 
ison to the previous observations with ASCA (Takano et 
al. 1994), BeppoSAX (Parmar et al. 2001), Chandra (La 
Parola et al. 2003) and XMM-Newton (e.g., Weng et al. 
2009), the signal statistics of the XIS spectrum, with a 
total FI and BI signal counts of 1.21 x 10 5 and 7.47 x 10 4 
respectively, seem to be considerably higher. The source 
was found to exhibit a featureless X-ray spectrum, without 
any clear sign of emission or absorption structure, except 
for instrumental and interstellar ones. 



For the purpose of quantifying the spectral properties 
of M33 X-8, we calculated the response matrix functions 
for the individual CCD chips with the HEADAS tool 
xisrmfgen, while the auxiliary response files were gen- 
erated by the tool xissimarf gen (Ishisaki et al. 2007), 
assuming a point source at the sky position of M33 X-8. 
The xissimarf gen tool precisely takes into account the 
reduction of the effective area due to contaminants on the 
optical blocking filters of the XIS, by referring to the ap- 
propriate version of the CALDB (sec §2). We performed 
the spectral fitting using the version 12.6.0q of XSPEC. 

As a first step, we tried to reproduce the observed XIS 
spectrum from M33 X-8 with an MCD model or a sim- 
ple PL one. A Galactic absorption was adopted with a 
free hydrogen column density. These models arc widely 
used as a standard tool for analyzing the X-ray spectra of 
Galactic black hole binaries and ULXs. As summarized 
in table 1, both models were found to be unacceptable 
( x 2 /dof = 3032.7/1373 and 2429.7/1373 for the MCD and 
PL models, respectively). The residual spectra with the 
MCD and PL models, shown in panels (b) and (c) of figure 
3 respectively, suggest that the observed X-ray spectrum 
of the source seems less convex than the MCD model, 
while it is more convex than the simple PL model. 

As a more sophisticated description for the X-ray spec- 
trum from an accretion disk, the BHSPEC model (Davis 
et al. 2005) 1 was fitted to the XIS data. The model 
assumes a fully rclativistic accretion disk around a spin- 
ning black hole, and self-consistently considers the vertical 
structure and radiative transfer within the disk. The disk 
spectrum was simulated for two representative values of 
the viscosity parameter, a = 0.1 and 0.01. There are fol- 
lowing four free parameters in the model; the black hole 
mass M, disk luminosity Ldisk, disk inclination angle i 

This model is available in the form com- 
patible with XSPEC from the following web 
http: //heasarc .gsf c .nasa. gov/docs/xanadu/xspec/models/bhspec 
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Fig. 3. Background-subtracted XIS spectrum of M33 X-8, shown without removing the instrumental response. Panels (a) and (e) 
compares the XIS data with the best-fit MCD+PL and extended MCD models, respectively. The residuals are displayed for (b) 
the best-fit MCD, (c) PL, (d) MCD+PL, (f) MCD+THCOMP, (g) extended MCD, and (h) slim disk (Kawaguchi 2003) models, 
respectively. 



(the face-on configuration corresponds to i = 0), and di- 
mensionless spin parameter a. Table 2 tabulates the resul- 
tant parameters. Although the fit was slightly improved 
by the BHSPEC model for a = 0.1 (x 2 /dof = 2118.6/1371) 
in comparison with the MCD or PL model, it was still in- 
adequate for interpreting the XIS spectrum of M33 X-8. 

Next, we examined a combination of an MCD com- 
ponent and a PL one (MCD+PL), both of which are 
subjected to a common free absorption. As shown in 
panels (a) and (d), a reasonable fit was derived by the 
model (x 2 /dof = 1535.8/1371) with the spectral param- 
eters listed in table 1. We obtained the absorption- 
inclusive observed X-ray flux as 1.57 x 10~ n ergs s _1 
cm" 2 in 0.5 - 10 keV. The absorption column density, 
A H = (1.80 ±0.13) x 10 21 cm -2 , became higher by a 
factor of 1.7 than the Galactic value (Ah = 1.1 x 10 21 
cm~ 2 ; Kalberla et al. 2005). By extrapolating the best- 
fit MCD component to the 0.001 - 100 keV range, we 
estimated the absorption-corrected bolometric disk lu- 
minosity as Ldisk = 2.9 x 10 38 (cosi) _1 ergs s _1 The in- 
nermost disk temperature of the MCD component was 
determined as T- m = 1.38^oog keV. The apparent inner- 
most disk radius is evaluated as r- m = 24.6+ 2 'g(cosi)~ a5 
km, based on the equation Ldisk = 47roT 2 n X^ (Mitsuda ct 
al. 1984; Makishima et al. 1986) where a is the Stefan- 



Boltzmann constant. By taking into account a factor £ to 
correct the inner disk boundary condition and a spectral 
hardening factor k representing the ratio of the color tem- 
perature to the effective one, r; n is converted to the true 
innermost disk radius as R la = t;K 2 r ln (Makishima et al. 
2000). According to the standard manner for the MCD 
model, we adopted £ = 0.412 (Kubota ct al. 1998) and 
k = 1.7 (Shimura & Takahara 1995). Thus, the true in- 
ner radius was evaluated as i?i n = 29.21?. '^(cosi) -0 5 km. 
We measured the photon index of the PL component as 
L = 2.301q;q7- The best-fit model indicates a dominance of 
the PL component over the MCD one, except for a narrow 
range of 3 - 5 keV. 

The PL component from the Galactic black hole bi- 
naries is thought to arise via Comptonization of the 
disk photons in the hot corona around the black hole 
and accretion disk. This motivated us to re-analyze 
the XIS spectrum by replacing the PL component with 
a thermal Comptonization continuum described by the 
nthcomp model in XSPEC (Zdziarski et al. 1999; here- 
after THCOMP). The seed photons for the THCOMP 
component were assumed to be supplied by the MCD one. 
As shown in panel (f) of figure 3, this MCD+THCOMP 
model became acceptable (x 2 /dof = 1549.6/1370), yield- 
ing the best- fit parameters tabulated in table 1. The 
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Table 2. Results from the BHSPEC or BHSPEC+THCOMP fitting 



Model 


BHSPEC ' 




BHSPEC+THCOMP 5 


N H (10 21 cm' 2 ) 


0.41 


0.25 


1.26tS-iU 


M (Mq) 


8.51 


13.0 


25.7 +2! ° 


cosi 


> 0.99 


0.67 


23 +ao7 

u - zo -0.03 


at 


> 0.79 


0.78 


a > 0.77 


i-disk (^e) 


0.80 


0.65 


40+ 14 


at 


0.1 


0.01 


0.1 


r 






9 91 +0.02 
^•■ £,L -0.11 


T sood (keV) 






n 91+O.O5 


T c (kcV) 






>4.5 



x7dof 2118.6/1371 2448.0/1371 1505.0/1367 



* All the errors are omitted, 
t Dimensionlcss spin parameter, 
t Fixed parameter in the model. 

§ The best-fit parameters are shown only for a = 0.1. 



model resulted in a full dominance of the THCOMP com- 
ponent in the XIS energy range, and the upper limit on the 
absorption-inclusive 0.5 - 10 keV flux of the MCD com- 
ponent was evaluated as Fx = 1.4 x 10 -13 ergs cm' 2 s" 1 
which corresponds to only ~ 1 % of the total X-ray flux. 
The disk temperature of the MCD component (i.e., the 
seed photon temperature) was found to be significantly 
lower than that in the MCD+PL model as T in = 0.36±o °4 
keV. We derived the asymptotic photon index and coro- 
nal electron temperature of the THCOMP component as 
r = 1.91+8$ and T c = 1.62±g;j$ keV, respectively. On the 
basis of the equation 



Y ' (Te/511 keV)[(r + 0.5) 2 -2.25] 

(Sunyacv & Titarchuk 1980), the combination of these 
two quantities indicates an optically thick corona with an 
optical depth of r ~ 15. 

We also considered a Comptonizing corona surround- 
ing an accretion disk around a spinning black hole 
by a sum of the BHSPEC and THCOMP components 
(BHSPEC+THCOMP). Because the BHSPEC model 
is not written in terms of the disk temperature, it 
is unable to constrain the seed photon temperature, 
Tseed, °f the THCOMP component to be consistent 
with the peak of the BHSPEC spectrum, through the 
BHSPEC+THCOMP fitting procedure alone. Actually, 
although a statistically reasonable fit was derived by the 
BHSPEC+THCOMP model with T sccd left free (see ta- 
ble 2), this model got settled to a physically unacceptable 
situation where the spectral peak of the BHSPEC com- 
ponent (> 4.5 kcV) highly surpasses the seed photon tem- 
perature (Tsccd =0.21^0 04 kcV). In addition, we failed to 
reconcile the seed photon temperature with the BHSPEC 
peak energy, in spite of a survey on the seed photon tem- 
perature in the range of T sco d = 0.3-2 keV. Therefore, in 
the present paper, we have given up to interpret the X-ray 
spectrum of M33 X-8 by Comptonization of X-ray photons 
from an accretion disk accompanied with a spinning black 



hole. 

Finally, in order to investigate the possibility of the slim 
accretion disk, which is regarded to be important at a high 
mass-accretion rate, we adopted a simple extension of the 
MCD model, where the radial temperature distribution on 
the disk is described as r~ p with the index p being a posi- 
tive free parameter (Mineshige et al. 1994). Although the 
model has been widely referred to as a p-free disk model, 
we regard this nomenclature as somewhat confusing. We, 
then, call it as the extended MCD model in the present 
paper. The model with p = 0.75 is identical to the simple 
MCD model, in which a standard accretion disk (Shakura 
& Sunyacv 1973) is adopted, while that with a smaller 
value of 0.5 < p < 0.75, is considered to approximate phe- 
nomenologically the X-ray spectra from a slim disk (e.g., 
Watarai et al. 2000). 

Panels (e) and (g) of figure 3 clearly visualizes that 
the extended MCD model reproduced the observed XIS 
spectrum of M33 X-8 well with a goodness of fit, 
X 2 /dof = 1541.8/1372, similar to that of the MCD+PL 
or MCD+THCOMP models. A rather flat tempera- 
ture profile with p = 0.5351^^5 was indicated by the 
model. We estimated the bolometric disk luminosity as 
Ldisk = 1-36 x 10 39 (cosi) _1 ergs s~ x . The innermost disk 
temperature, T- m = 2.00lo'o5 kcV, was found to be higher 
than that in the MCD+PL model. In the similar man- 
ner as for the MCD+PL model, we obtained the apparent 
innermost disk radius as r; n = 25.8^2 '^(cosi) -0 - 5 km. In 
the case of the slim disk, the spectral hardening factor is 
suggested to become larger than that in the standard ac- 
cretion disk (k = 1.7; see above), due to various effects in- 
cluding enhanced electron scattering and Comptonization 
(e.g., Watarai & Mineshige 2003; Kawaguchi 2003). Here, 
we adopt tv — 3 clS cL typical value for the slim disk, which 
was required to reconcile the observed high disk temper- 
ature of the Galactic microquasar GRS 1915+105 with 
the prediction from the numerical simulation on the disk 
at a high accretion rate (Watarai & Mineshige 2003). 
In addition, the boundary condition is corrected by as- 
suming £ = 0.353, which was derived by taking into ac- 
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count the transonic flow in the pseudo-Newtonian poten- 
tial (Vierdayanti et al. 2008). These yielded the "true" 
disk radius as R in = 81.9±^(cosi)"°- 5 (£/0.353)(K/3) 2 km. 

4. Discussion 

4-1- Summary of the Suzaku observation 

In the Suzaku observation conducted on 2011 January 
11 - 13 with an effective exposure of 91.0 ks, we success- 
fully acquired the X-ray spectrum of the nearest ULX, 
M33 X-8, with the highest signal statistics in the 0.4 - 
10 keV range. During the observation, the object exhib- 
ited no significant change in its spectral hardness, with 
only a weak intensity variation with a standard devia- 
tion of ~ 5%. Without being corrected for the absorp- 
tion, the time-averaged 0.5 - 10 keV flux of the source 
was measured as 1.57 x 10~ n ergs s _1 cm -2 (evaluated 
from the best-fit MCD+PL model). Re-evaluated from 
the previous results (Takano ct al. 1994; Parmar et al. 
2001; La Parola et al. 2003; Foschini et al. 2004; Weng et 
al. 2009; Middleton et al. 2011), its 0.5 - 10 keV flux was 
found to stay within a relatively narrow range of (1.4 - 
1.9) x 10 -11 ergs s _1 cm -2 over nearly 16 years 2 . Thus, 
the source X-ray flux during the Suzaku exposure is lo- 
cated in the middle of its distribution in the previous ob- 
servations. 

Figure 3 illustrates that the time-averaged XIS spec- 
trum of the source was described by the MCD+PL, 
MCD+THCOMP, or extended MCD one, with the best- 
fit parameters summarized in table 1, while the single 
PL, MCD or BHSPEC model was found to be unaccept- 
able. Assuming the best-fit extended MCD model, the 
XIS spectrum was de-convolved into the spectral energy 
distribution in the vF v form, as shown in figure 4. These 
results are utilized to infer the nature of M33 X-8 in the 
following. 

4-2. Comptonization scenario 

When we adopt the MCD+PL model, both the disk 
temperature and the PL photon index of M33 X-8, T- m = 



Table 3. Physical parameters evaluated by the Slim disk 
model (Kawaguchi 2003) 



1.38±°,;^ keV and T 



2.30± ;^, are found to fall within 



the range derived in the previous results (e.g., Takano et 
al. 1994; Parmar et al. 2001; La Parola et al. 2003; Weng ct 
al. 2009). Especially, the steep photon index appears to be 
close to those of Galactic black hole binaries in the very 
high state (r > 2.4; e.g., Miyamoto et al. 1991; Kubota 
& Done 2004). Actually, the very-high-state interpreta- 
tion was proposed for the source by Foschini et al. (2004), 
based on the XMM-Newton spectra obtained in a few oc- 
casions. However, we regard that this interpretation en- 
counters several difficulties. 

The first obstacle for the very-high-state scenario comes 



2 Here, we neglected several XMM-Newton observations (ObsID 
= 0102642001, 0141980501, and 0141980601), for which we 
found inconsistency of more than 20 % in the flux determination 
between the separate results (Foschini et al. 2004; Weng et al. 
2009; Middleton et al. 2011). Even if these results are taken into 
account, the observed 0.5 - 10 keV flux varied by only a factor 
of ~ 2 in the range of (1.1 - 2.4) X 10 -11 ergs s _1 cm~ 2 . 



Parameter 



Value 



A H (10 21 cm" 2 ) 
M (M ) 
M (L E /c 2 ) 
a 

Lx (ergs s" 1 ) * 
£disk (ergs s" 1 ) f 



1.23] 



RTUB 

'-0.05 

13.11°;? 
10.5t°; 5 7 
0.21 + 0.05 
0.68 x 10 39 
0.94 x 10 39 



yVdof 1516.7/1372 

The model assumes a tacc-on disk configuration 
(i = 0). 

* Absorption-corrected 0.5 - 10 keV luminosity, 
t Absorption-corrected bolomctric disk luminosity 
evaluated in 0.001 - 100 keV. 



from the high disk temperature of M33 X-8, which sig- 
nificantly exceeds that of Galactic black-hole binaries 
in the strongly Comptonized very high state (T; n <C 1 
keV), termed intermediate hard state. In the case of the 
Galactic black-hole binaries, such a low disk temperature 
in the very high state is explained by a possible disk trun- 
cation at a radius significantly larger than the last sta- 
ble circular orbit (Kubota & Done 2004). Secondly, the 
PL dominance at energies below the disk temperature T- m 
seems quite unphysical. The PL component in the very 
high state is naturally predicted to exhibit a sharp low- 
energy cut-off below the seed photon energy, because it 
is basically attributed to Comptonization of disk photons 
within the hot corona, surrounding the black hole and 
accretion disk. 

We tried to reconcile these inconsistencies, by replac- 
ing the PL component with the Comptonization contin- 
uum described by the THCOMP mode, where the MCD 
emission was assumed to provide the seed photon source. 
The MCD+THCOMP model has eventually arrived at a 
different solution that the disk photons are Comptonized 
by a low-temperature (T c = 1.621q o? keV) optically-thick 
(r ~ 15) material, which is clearly deviating from the very 
high state (T c > 20 keV and r < 2; Kubota & Done 2004). 
Recent numerical simulations anticipate that extremely 
high radiation pressure at a highly super-Eddington mass 
accretion rate (M ^> Le/c 2 ) induces a strong outflow 
from the disk (Ohsuga et al. 2009), which serves a cool 
(T e <C 10 keV) optically thick (r>l) Comptonizing corona 
(Kawashima et al. 2009). Such a spectral features are 
gradually suggested from several high-luminosity ULXs 
(e.g., Gladstone et al. 2009; Vierdayanti et al. 2010). In 
addition, Middleton et al. (2011) applied a similar idea 
to the spectral variation of M33 X-8. However, the high- 
quality XIS spectrum from the source requires the ab- 
sorption column density for the MCD+THCOMP model 
(A H = 0.85l°;i4 x 10 21 cm" 2 ) to be lower than the Galactic 
value (A H = 1.1 x 10 21 cm" 2 ; Kalberla et al. 2005). We 
regard such a situation invalid, considering the absorp- 
tion within the M33 galaxy and intrinsic to the source. 
Therefore, we have concluded that the Comptonization is 



8 



N. Isobe, et al. 



Spectral Enedgy Distribution ofM33 X-8 




Fig. 4. Absorption-uncorrccted spectral energy distribution 
of M33 X-8. The XIS data were de-convolved with the best-fit 
extended MCD model. 

not a realistic scenario for M33 X-8. 
4-3. Slim-disk interpretation 

From observational (e.g., Kubota & Makishima 2004) 
and theoretical (e.g., Watarai et al. 2000) studies of 
Galactic black hole binaries and their accretion disk in the 
decade, the extended MCD model is recognized as a good 
approximation to the X-ray spectrum from a slim accre- 
tion disk, which is thought to be realized at a high X-ray 
luminosity. Actually, the model was successfully utilized 
to interpret the X-ray spectra of Galactic black hole bina- 
ries, at the highest end of its X-ray luminosity or in the 
sub-/trans-Eddington regime, as the slim disk state (e.g., 
XTE J1550-564; Kubota & Makishima 2004). These were 
subsequently followed by successful applications to the X- 
ray spectral properties of several ULXs (e.g., Tsunoda et 
al. 2006; Mizuno et al. 2007; Isobe et al. 2009). 

We demonstrated that the XIS spectrum of M33 X-8 
was successfully described by the extended MCD model. 
The innermost disk temperature and radial temperature 



gradient, T in = 2.001^ kcV and P = °- 535 -o:oo5 respec- 
tively, were both found to agree fairly with those previ- 
ously inferred from the object (T; n = 1.4 - 2.5 keV and 
p = 0.50 - 0.56; Weng et al. 2009). Especially, the flat ra- 
dial temperature distribution within the disk agrees with 
the prediction from the slim-disk theory (Watarai et al. 
2000). Therefore, the slim-disk picture seems convincing 
for the source. 

Different from Galactic black hole binaries in the classi- 
cal high/soft state, which have been observationally recog- 
nized to harbor a standard accretion disk extending down 
to the innermost stable circular orbit at 3i?s, a mass 
estimation from the assumption of i?j n = 3i?s dose not 
hold for a black hole with a slim disk, where i?s is the 
Schwarzschild radius. This is because theoretical simula- 
tions on the slim disk widely predict that emission from 
the accretion flow inside the innermost stable circular or- 
bit becomes not negligible (Watarai et al. 2000). In rela- 



tion, apparent dependence of the innermost disk radius 
(and hence the luminosity) upon the disk temperature 
is suggested from the numerical simulation (Watarai et 
al. 2000) as i? in ~ T^ 1 (corresponding to i disk ~ T£). 
Similar properties were actually observed from several 
ULXs within the framework of the standard MCD+PL 
spectral modeling (e.g., Mizuno et al. 2001). 

In order to cope with this problem, Vierdayanti et al. 
(2008) proposed a simple method. They fitted the nu- 
merically simulated X-ray spectra of a slim disk to the 
extended MCD model. They found from some represen- 
tative values of black hole mass, M = (10 100)M Q , and 
accretion rate, M = (1 - 1000)Le/c 2 , that a correction 
factor needs to be applied to convert the apparent mass 
Mx derived from the hypothesis of R m = 3i?s to the true 
one as M/Mx = 1.2 - 1.6. Here, we follow the technique to 
estimate the black hole mass of M33 X-8, assuming that 
the object accompanies a slim accretion disk. 

The innermost disk radius of M33 X-8, i?i„ = 
8L9tg5(cosi)~ - 5 km (here and hereafter dependence on 
k and £ is neglected for clarity), measured from the XIS 
spectrum with the extended MCD fitting gives the ap- 
parent black hole mass as Mx = (9.2 ± O.7)(cosi) _O5 M0. 
Adopting the conversion factor of M/Mx ~ 1.2 (for M = 
1OM ; Vierdayanti et al. 2008), the true mass of M33 X- 
8 is estimated as M ~ ll(cosi)~ a5 M Q . Since the mass 
corresponds to the Eddington luminosity of Le ~ 1.7 x 
10 39 (cosi) -0 5 ergs s _1 , the source is indicated to shine at 
the sub-Eddington luminosity of Lx/Le ~ 0.4(cosi)~ - 5 
or idisk/^E ~ 0.8(cos i)~°- 5 . Considering the observa- 
tional fact that the Galactic black hole binaries are known 
to exhibit the slim-disk-likc X-ray properties typically at 
the Eddington ratio of Lx/Le = 0.3 - 1 (e.g., Kubota & 
Makishima 2004; Kubota et al. 2001a), the slim-disk in- 
terpretation is thought to become self-consistent for M33 
X-8. 

A number of authors (e.g., Okajima et al. 2006; 
Vierdayanti et al. 2006; Yoshida et al. 2010) utilized the 
X-ray model spectrum from a slim accretion disk gener- 
ated from the numerical simulation by Kawaguchi (2003) 
3 to claim a super-Eddington interpretation for several 
ULXs, including M33 X-8 (Foschini et al. 2006). These 
have driven us to investigate the XIS spectrum with this 
model. The model takes into account the effects of elec- 
tron scattering and general relativity in detail (with the 
parameter Model ID set at 7). We have to note that only 
the face-on disk geometry (i = 0) is supposed in the model. 
By leaving free M, M and the viscosity parameter a, the 
model yielded a reasonable fit (x 2 /dof = 1516.7/1372) as 
shown in the panel (h) of figure 3. The physical quan- 
tities derived from the slim disk model are tabulated in 
table 3. The black-hole mass, M = 13.lt°;?M Q , was found 
to nearly agree with the result from the extended MCD 
model. The model requires the viscosity parameter of 
a = 0.21 ± 0.05, which is not deviating far from that usu- 



This model is available in the form com- 
patible with XSPEC from the following web 
http: //heas arc .gsfc .nasa. gov/docs/xanadu/xspec/models/slimdisk. html. 



No. ] 

ally adopted for ULXs (a = 0.1; e.g., Vierdayanti et al. 
2006; Yoshida et al. 2010). Foschini et al. (2006) reported 
to prefer a lower mass of M < IOMq, by applying the same 
model to the XMM-Newton spectra of M33 X-8 obtained 
in the three occasions. However, we regard that their re- 
sult is possibly an artifact due to a rather low value of 
the disk viscosity (a < 0.01) adopted in their analysis. It 
is important to note that the luminosity of the ULX is 
indicated to be significantly less than the Eddington limit 
as Lx/Le ~ 0.3 or Ldisk/^E ~ 0.5, even though the source 
is suggested to be accreting at a super Eddington rate as 
M=l0.5±° - 5 7 L E /(?. 

We regard these results as basically supportive of the 
slim disk idea with the sub-Eddington luminosity for M33 
X-8. However, we have to be careful about a few the- 
oretical and observational uncertainties. The black hole 
mass derived from the extended MCD model is subject 
to the ambiguity in the spectral hardening factor k. The 
numerical model on the slim accretion disk (Kawaguchi 
2003) is recommended to be calibrated by the Galactic 
black hole binaries with a known mass, although they are 
known to exhibit rarely the slim disk properties. In ad- 
dition, it is expected from recent numerical studies (e.g., 
Ohsuga et al. 2009; Kawashima et al. 2009) that the slim 
disk tends to accompany a strong outflow. The X-ray 
spectrum from the slim disk is possible to be contami- 
nated from Comptonizcd X-ray photons produced in the 
outflow. 

5. Implication to the nature of the ULXs 

By making most of the Suzaku XIS spectrum with the 
highest signal quality ever recorded, we propose the slim 
disk scenario for the nearest ULX, M33 X-8. In addition, 
both the empirical extended MCD model and the com- 
plicated numerical simulation on the slim accretion disk 
have guided us to the conclusion that the source hosts a 
black hole with a mass of M > 10M S (for i = 0). If the 
disk viewing angle is taken into consideration, the mass 
will become larger as oc (cosi) -0 ' 5 in the simplest manner 
neglecting the effects of the general relativity. One of the 
most remarkable points from this investigation is that the 
source is found to be radiating at the sub-Eddington lumi- 
nosity, in spite of its super-Eddington mass accretion rate. 
A similar observational support was revealed by Isobe et 
al. (2009), to the sub- or trans-Eddington luminosity from 
the slim accretion disk by analyzing the Suzaku data for 
the ULX, Source 3 in the nearby galaxy NGC 2403. 

A sub/trans-Eddington luminosity is not surprising 
from the theoretical point of view, even in the case of the 
super critical accretion flow (M S> Le/c 2 ), since the lumi- 
nosity tends to saturate due to strong advection and pho- 
ton trapping. Actually, a disk luminosity of Ldisk/LE ~ 3 
at most was anticipated from the numerical simulation on 
the slim disk at the significantly super critical accretion 
rate as M = 1000L E /c 2 (e.g., Vierdayanti et al. 2008). In 
summary, we have asserted that it is not necessary to in- 
voke the highly super-Eddington luminosity, in order to 
understand the ULX behavior. By extrapolating this con- 
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sidcration to the ULXs more luminous than M33 X-8, it 
is natural to conclude that there are intermediate mass 
black holes with a mass of M > IOOMq among those with 
Ldisk> 10 40 ergs s" 1 . 

The bandpass of most operating X-ray instruments lim- 
ited below 10 keV frequently prevented us to disentan- 
gle clearly the degeneracy between the candidate spectral 
models (e.g., the slim disk or cool Comptonization model). 
Actually, the capability of the hard X-ray spectrum was 
certified by the recent Suzaku HXD detection of the hard 
X-rays up to ~ 20 keV from the ULX, M82 X-l (Miyawaki 
et al. 2009), since the result have succeeded in specifying 
the source spectrum as a Comptonized emission from a 
black hole with a mass of M = (100 - 200)M Q . Moreover, 
even in the case that a source is confirmed to host a slim 
accretion disk, a hard band X-ray spectrum above 10 keV 
is expected to be useful to eliminate the possible uncer- 
tainty due to the contamination from the strong outflow 
(see §4.3). Therefore, we expect that a wide-band X-ray 
spectrum in the range of 0.3 - 80 keV to be derived with 
the next-generation Japanese X-ray observatory ASTRO- 
H (Takahashi et al. 2010) will be inevitably play a crucial 
role to settle the nature of the ULXs. 

We thank the anonymous referee for her/his supportive 
guidance to improve remarkably the present paper. We 
are grateful to all the members of the Suzaku team, for 
the successful operation and calibration. This research is 
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Scientists (B) 22740120 (N. L). 
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